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Abstract 

We give an exact information-theory treatment of the n-dimensional black- 
body radiation process in a non-extensive scenario. We develop a q gener- 
alization of the laws of i) Stefan Boltzmann, ii) Planck, and iii) Wien, and 
show that conventional, canonical results are obtained at temperatures above 
IK. Classical relationships between radiation, pressure, and internal energy 
are recovered (independently of the q- value). Analyzing the particles' density 
for q ~ 1 we see that the non-extensive parameter q introduces a fictitious 
chemical potential. We apply our results to experimental data on the cosmic 
microwave background and reproduce it with acceptable accuracy for different 
temperatures (each one associated to a particular q value). 

PACS: 05.30.-d, 05.30.Jp 

KEYWORDS: g-Thermostatistics, black-body radiation. 



martinez @ venus .fisica. unlp . edu . ar 
pennini@venus. fisica. unlp.edu. ar 
plastino@venus .fisica. unlp .edu. ar 
tessonec@venus.flsica.unlp. edu. ar 



1 



I. INTRODUCTION 



Black-body radiation studies constitute a milestone in the history of Physics. Planck's 
law is satisfactorily accounted for by recourse to Bose-Einstein statistics and has been ex- 
perimentally re-confirmed over and over for almost a century. However, in the last decade 
small deviations from this law have been detected in the cosmic microwave radiation M . In 
0, its authors advance the hypothesis that these deviations could have arisen at the time 
of the matter-radiation de-coupling, and associate them to a non-extensive statistics envi- 
ronment. They argue in favor of such an hypothesis on the basis of the close relation that 
exists between long-range interactions and non-extensive scenarios Hf§]. The concomitant 
non-extensive thermostatistical treatment f3j-|T2"H is by now recognized as a new paradigm for 
statistical mechanical considerations. It revolves around the concept of Tsallis' information 
measure S, |, a generalization of the logarithmic Shannon's one. S q is parameterized by a 
real index q and becomes identical to Shannon's measure when q = 1. 

The study reported in employs the so-called Curado- Tsallis un-normalized expectation 
values in the limit q — > 1. Nowadays, it is believed that the Curado-Tsallis (CT) frame- 
work has been superseded by the so-called normalized or Tsallis-Mendes-Plastino (TMP) 
II] , II] , which seems to exhibit important advantages [[/J. This normalized treatment, 



one 



in turn, has been considerably improved by the so-called "Optimal Lagrange Multipliers" 
(OLM) approach |15| . It is, then, natural to revisit the problem from such a new viewpoint. 

The OLM treatment can be recommended in view of the findings of |l6 Hl8[l regarding 
the particular nature of the Lagrange multiplier associated to the temperature. Within an 
OLM context, non-extensivity is restricted just to the entropy. The internal energy remains 
extensive and the Lagrange multipliers conserve their traditional intensive character (al- 
lowing one to identify them with their thermodynamic counterpart). Moreover, the OLM 
approach unifies Tsallis and Renyi variational formalisms under a g-thermostatistics um- 
brella. The solutions to the concomitant OLM-Tsallis variational problem are also solutions 



for its Renyi-counterpart |19|. Since Renyi's entropy is extensive, one easily understands 
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thereby the OLM formalism's success in reproducing classical thermodynamic results in a 



simple and clean fashion [20.21 



We shall study here, from an OLM viewpoint, black-body radiation in equilibrium within 
an enclosure of volume V, with the goal of ascertaining the possible g-dependence of i) the 
Planck spectrum, ii) Stefan-Boltzmann's law, and iii) Wien's one. We will compare the 
ensuing results with those found in the literature Pj22|-p4| (all of them under Curado-Tsallis 
treatment) and with experimental data [l|]. Due to the fact that most previous non-extensive 
treatments of quantal gases employ the so-called Factorization Approach (FA), and that 
this approximation underlies some of the preceding treatments of the black-body radiation 
problem, we will translate the FA into the language of a normalized-OLM Factorization 
Approach (OLM- FA). This in turn provides one with a g-normalized generalization of the 
particle's density expression for quantum gases. 

The paper is organized as follows: In Sect. we present a brief OLM primer. In Sect. 
PTI| we obtain the OLM partition function. Sect. [IV] is devoted to the Stefan-Boltzmann 
law. We obtain the exact expression for the internal energy ( |1V A| ) and find that the Stefan- 
Bolzmann law holds for all temperatures except a small interval (T e [lO -2 ^, IK]). We also 
revisit pertinent antecedents in the literature of the non-extensive approach to the black- 
body problem (|I V B[ ) . In Sect. [V| we perform an exact calculation of the energy density. 
Planck's law (|V A ) seems to be valid everywhere, with just slight variations in the shape 
of the associated curves. Following ||, the formalism is applied to the cosmic microwave 
radiation ( |VQ) in order to search for deviations from Planck's law. Using the maximua 
of the energy density curves, we search for the g-generalization of Wien's law ( |V B|) and 
find good agreement, save for a small temperature interval (T £ [10~ 4 K,10~ 2 K]). Some 
conclusions are drawn in Sect. |VJ. In Appendix [A| we review the FA, introduce our new 
OLM-FA technique and apply it to the black-body radiation problem. 
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II. MAIN RESULTS OF THE OLM FORMALISM 



For a most general quantal treatment, in a basis-independent way, consideration is re- 
quired of the density operator p that maximizes Tsallis' entropy 

S q 1 - Tr(pi) 



k 



q-l 



subject to the M generalized expectation values (Ojj , where Oj (J — 1, 
relevant observables. 



(1) 

M) denote the 



Tsallis' normalized probability distribution [n| is obtained by following the well known 
MaxEnt route |25 |. Instead of performing the variational treatment of Tsallis-Mendes- 
Plastino (TMP) ||13|| , we will pursue the alternative path developed in [15|. One maximizes 
Tsallis' generalized entropy given by Eq. (U) |||6]|| subject to the constraints |5|JT5| 



Tr{p) = 1 



Tr 



p q [o 3 -(o. 



0. 



(2) 
(3) 



whose generalized expectation values 



Tr(p q Oj) 
Tr(pi) 



(4) 



are (assumedly) a priori known. In contrast with the TMP-instance [ 13[| , the resulting 
density operator 

fl/(l-q) 



P 



(5) 



is not self referential [|T^|. In Eq. fl5|), {\j} (J = 1, ...,M) stands for the so-called "optimal 
Lagrange multipliers' set". The quantity f q is known as the configurational characteristic 
and takes the form 



M , v 
/, = l-(l-ff)E^(Oi-(Oi)J, (6) 
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if its argument is positive, while otherwise f q = (cut-off condition |T3|]). Of course, Z q is 
the partition function 



z q = T r m-i\ 



(7) 



It is shown in Ref. |T3J that, as a consequence of the normalization condition, one has 



K q = Trf*^ = Z { 



gi 



(8) 



which allows one to write Tsallis's entropy as 



S q = k \n q Z q , 



(9) 



with ln q x = (1 — x 1 ~ q )/(q — 1). These results coincide with those of TMP ||13|| . Using Eq. 
(0), the connection between the OLM's set and the TMP's Lagrange multipliers set can be 
written as [[OJ 



X TMP = 



(10) 



The TMP's Lagrange multipliers are not intensive quantities fiqjlq| . They do not, as a 
consequence, have a simple physical interpretation. To the contrary, the OLM multipliers 
are intensive |16| , [l8|| , a fact that can be easily explained noticing that they are the natural 
Lagrange multipliers of a Renyi's variational approach |19| . 

If we define now 



(11) 



we are straightforwardly led to [16 







(lnZ ? ) =\j, 



d 



dXj ■ > ' ■/ 



O, 



(12) 
(13) 



Equations ( |T2"D and ( [13] ) are fundamental information theory (IT) relations for formulat- 
ing Jaynes' version of statistical mechanics [F25J. Due to Eqs. (|S[) and (]TI1|), the IT relation 
(|l2l) leads straightforwardly to the well known expression [|L| 







d O, 



~q \ _ yTMP 



(14) 
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III. PARTITION FUNCTION AND RADIATION PRESSURE 



We will introduce now the exact (OLM) black-body radiation treatment. The chemical 
potential is taken, of course, equal to zero (Grand Canonical ensemble with \x = 0) and, in 
looking for the equilibrium properties we consider that the appropriate thermodynamical 
variables are, as customary, the volume V, and the temperature T Consider first the 
standard situation q = 1. 

The Hamiltonian of the electromagnetic field, in which there are n^ e photons of momen- 
tum k and polarization e, is given by 

H = J2huh kje , (15) 

k,e 

where the frequency is u = c |k| and Uk, t = 0, 1, 2, . . ., with no restrictions on {^k,e} • 

For a macroscopic volume V the density of states g in an n-dimensional space is g n {uj) = 
A n Lj n -\ with 

2r n V 

n (47rc 2 )«/ 2 r(n/2)' 1 ' 

where r n = n — 1 is the number of linearly-independent polarizations. 
The partition function 



Z x = Tr (e- 1 ™) (17) 

can be written as 

Z x = exp jy dujg n (uj) In [1 - exp J = e ? ", (18) 



where 



T A 

± n- rl n 



(19) 



and 



I n = - r dxx"- 1 In (l - e~ x ) = r(n)C(n + 1), (20) 

J 



with ( standing for the Riemann function and T for the Gamma function. 

The OLM-Tsallis generalized configurational characteristic will be (Cf. Eq. (||)) 

f q = l-(l-q)P(H-U q ), (21) 

where H is the Hamiltonian given by Eq. (|T5|), U q is the mean energy introduced in Eq. 
(1), and (3 = 1/kT. 

With the aim of calculating Z q as defined by Eq. (0), we follow the steps of Ref. f22 



and use the integral (Gamma) definition given by the relation |29 



7 dt m^ for 6 > 
27r -oo (22) 

for 6 < 0, 

with Re(z) > 0, and —n/2 < arg(a + it) < ir/2. 

If we set b = f q (the cut-off condition is naturally fulfilled J!]]) and z = 1/(1 — q) + 1 
(Re(z) > 0, so that either q > 2 or q < 1), the generalized partition function adopts the 
appearance 

oo 

Z q (U q ) = J dtK q (t)Z x {p), (23) 

with Zi given by Eq. 



-oo 



K m _ r[(2-g)/(l-g)]exp(l + zt)e^ 

27T(l + Zt)( 2 ^)/(^) ' 

and 

/3=(l + zt)(l-g)/?. (25) 
In order to evaluate the integral in Eq. (|23|), we expand the exponential and obtain 

Z q {U q ) = — , jJ E %(1 - / Ei— > (26) 

27r ^0™ ! (l + ^t)T4+" m 



where £ n is given by Eq. (|19f) . 

Using again Eq. (p2|), with 6 = 1 + (1 — q)(3U q and z = (2 — q)/(l — q)+ nm, we arrive at 



oo 

Z q (U q ) = r [(2 - q)/(l - q)] [1 + (1 - q)PU q ] l/{1 ~ q) E t( 2 " ?)/(! " ?) + *H , ( 2 ?) 



m=0 

where 



1 77 m 



Notice that an additional, cut-off-like condition must be considered, namely, 

1 + (1 - q)PU q > 0, otherwise, Z q (U q ) = 0. 

A similar path can be followed in order to obtain the quantity lZ q (introduced in Eq. 
(§)), which will read 

oo 

K q {U q ) = r [1/(1 - q)} [1 + (1 - q)PU q } 9/{1 ~ 9) E Brn^ 1 [1/(1 ~q) + H , (29) 

m=0 

although in this case the allowed interval of q- values is reduced to < q < 1. Note that this 
permissible interval of g-values respects the new cut-off condition introduced above. 

Once we have Z q , the radiation pressure is easily obtained by applying Eq. (|l2f) to 
(PP,V) H,i.e., 

, p _l, 7 l EmmB m r- 1 ((2-q)/(l-q)+nm) 

Pr dV Zq V E m B m V~^ ((2 - g)/(l - g) + rnn) ' 1 J 



Due to the fact that we have set /i = 0, the quantity Z q of Eq. fl27|) does not depend on 
the mean number of particles iV ? . One has 

/* = j 9- 1 ^ln^ = 0. (31) 

IV. STEFAN-BOLZMANN'S LAW 

A. Exact OLM treatment 

As in the previous section, the generalized internal energy of the black-body radiation is 
obtained from Eq. (Q) by specializing the problem to the Grand Canonical ensemble with 
fi = 

8 



U q = K- l Tr {jf^H) , (32) 



with f q given by Eq. (pT|). 

Consider now the trace's content. Using (^) we can find 

,(l+it)[l+(l-q)f3U q ] 



, rn/fi — oil 7 e v+mi+{i- q )pu q \ _ 
tf 9 = V l / o / ^— z^ Tr ( Z M H ) > (33) 



where /3 is given by Eq. (p5|). The permissible g- interval is < g < 1, due to restrictions 
posed by the Gamma integral representation. Taking advantage of the fact that 



Tr 

we obtain 



J 3/3 (StT ' 



OO 



n / 1 \ £ n 00 £™ I 7 e (i+it)[i+(i-?)/Jtf<,] 



^ ^^/(l-gj (l-^S^Kl-ff^i, dt (l + iO^ +n(1+m)+1 ' (35) 
and, using again Eq. fl22|), we obtain the mean energy expression 



„t E 5 m r- 1 [l/(l - g) + ra(m + 1) + 1] 

U q = n [1 + (1 - ff) !I ^o • (36) 

(1 - g)"+^ E 5 TO r- 1 [l/(l - g) + nm] 



Notice that the Tsallis'cut-off condition 1 + (1 — q)(3U q > is always satisfied. The series 
in Eq. ( pop rapidly converges on account of the exponential-like factors B m (Cf. Eq. fl28|)), 
weighted by inverse Gamma functions. 

Lenzi et al. solved the g-black-body radiation problem in exact fashion in its Curado- 
Tsallis un- normalized version ||22|| . Since their resulting internal energy is not self referential 
(as it is in the TMP normalized instance) they were able to describe the asymptotic behavior 
for (3 in quite simple terms. Here we need to perform a more detailed analysis by considering 
different possibilities for the form that the product (3U q may take as (3 — > oo. For instance, 
if we assume that i) (3U q goes over to a constant or ii) it is not bounded, the limiting process 
leads to incoherencies, while if we assume that, in Eq. (^), (3U q — > when (3 — > oo, U q oc T 4 , 
as one has the right to expect. 



9 



As a consequence of the normalization condition given by Eq. ([3D, we know that TZ q = Z q 
(see Eq. (pi)). This relation allows us to look for alternative expressions for the relevant 



mean value. If we evaluate U q in terms of Z q we find 



E mB m T- 1 ((2 - g)/(l - g) + nm) 

^ = • (37) 

P E 5 m r- 1 ((2-g)/(l-g)+nm) 



m=0 



By inspection of Eq. (30) one realizes that the traditional relation between internal 



energy and pressure still holds here 



P=I£. (38) 



Eqs. ( p6| ) or (|3~T| ) are recursive expressions that give the OLM version of Stefan- 
Boltzmann's law, which, for the q = 1 case, reads 

U = ninjfi oc T n+1 . (39) 

The present equations are to be tackled numerically. Fig. [I] depicts U q as a function of T for 
different values of q and n = 3 in a log-log scale, where {7 9 has been evaluated from Eq. (p6|). 
It is seen that Stefan-Boltzmann's law is reproduced by our formalism for a wide range of 
T- values. Violations are detected just for some special T-ranges (that depend on q), within 
the interval lO" 2 ^ < T < IK. 

With the present formalism the Stefan-Bolzmann constant becomes a function of i) q 
(namely, u a q ) and ii) the relevant range of temperatures. Indeed, for temperatures below 
lO -2 -^ results will be markedly different from those obtained if we consider T > IK (see 
Fig. 0). It may be appreciated that, for these two Temperature ranges, the values of a q 
increase monotonically with q. 

Looking for insights into the meaning of the "violation range" lO -2 ^ < T < IK we 
considered the specific heat 

C. = f. (40, 
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The concomitant results are plotted in Fig. [3|. One sees that the violation range coincides 
with that of a constant specific heat, i.e., U q oc T. The specific heat curves' behavior is 
quite different from that associated with Gibbs' predictions, even for q values close to unity. 
By recourse to detailed numerical analysis one notes that the peak one observes is not a 
discontinuity but part of a smooth curve that reflects the typical behavior of a first excited 
energy level when that level lies too close to the ground state. The typical step-like form 
tends to disappear in the q — > 1 limit, just as if 1 — q were a new "degree of freedom" of the 
system. 

For more details on the transition between Gibbs' and Tsallis's statistics let us analyze 
the q — > 1 limit in Eq. ([36]). Using the fact that the pertinent series are highly convergent, 
we can try to recover analytically the above numerical results by keeping only the first order 
term of the pertinent series, 

U q « ^ [1 + (1 - q)/3U q ] n+1 , (41) 

where Xg — T[l/(1 — q)](l — g) _ ( n+1 'r -1 [l/ (1 — q) + n + 1]. For q — > 1, Xg — > 1 an d we see 
that Stefan-Boltzmann's law is recovered. 
Now, to first order in 1 — q, 



[1 + (1 - q)(3U q ] n+1 « 1 + (n + 1)(1 - q)pU< 



so that, rearranging terms, Eq. (fTI]) can be cast as 



q I - (I - q)( n + l)n Xq in 1 1 

The change of behavior we are interested in is better observed with reference to Eq. (|37|), 
by keeping only first order terms in the series expansion. We find 

rr „ n B 1 T- 1 [(2-q)/(l-q) + n] 

q ~ P^[(2-q)/Q.-q)) + B 1 r-i[(2-q)/(l-q)+nY 1 J 



where 



+ (44) 
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It is clear that, for j3U q such that the first term in the denominator of Eq. ( f43|) dominates, 
U q oc T n+1 . When the second term is dominant, instead, then U q oc T, in agreement with 
our numerical results. Eq. ( [12] ) displays a similar behavior, although its validity is restricted 
to the q — > 1 limit. This prevents the second term in the denominator from being dominant, 
a fact reflected in Fig. [I], where the linear dependence is seen to fade away. 



B. Comparison with previous non-extensive results 

The Stefan-Boltzmann's law was first discussed within the Tsallis' non-extensive frame- 
work in p|. In this paper the authors employ the so-called Curado-Tsallis un-normalized 
expectation values || and work in the limit q — > 1: a first order approximation in 1 — q 
is used for the partition function in order to study the cosmic microwaves' background (in 
Sect. |V C we will apply the results of the present work to the same data set). In a subse- 



quent effort, Lenzi et al. |22|] advanced an exact treatment for the same problem, also within 



the Curado-Tsallis framework. A third relevant work is that of Tirnakli et al. |24j], that 
compared the exact nonextensive treatment of the problem with the FA one, including the 
un-normalized asymptotic approach (AA) already introduced in ||. 

The different ensuing expressions for the pertinent a la Stefan-Bolzmann laws are given 



in Table |A 2j for n = 3. The traditional Gibbs expression can be found in the first row, while 
the second and third are, respectively, the Curado-Tsallis' and OLM (exact) results. The 
TMP solution can be read off the fourth row, where we have simply used Eq. ([H]) for (3. 
The last three rows are devoted to the q — > 1 limit. The first one contains the factorization 
approach result, the second is the our new, OLM factorization approach expression, and, 
finally, the last one gives the exact OLM treatment in this limit. 

Inspection of these expressions allows one to appreciate that the OLM-FA seems to 
inherit characteristics of both the FA and the OLM treatments for gwl. The OLM-FA 
numerator coincides with that of the FA approach, while the denominator resembles the 
OLM one (except for the sign). Indeed, the concomitant integrals are identical here to those 
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of the un- normalized treatment. 



V. PLANCK AND WIEN LAWS 



A. Planck's law: exact OLM treatment 



The generalized spectral energy distribution of black-body radiation u q is defined by the 
integral 

oo 

U q — J duJU q . (45) 


In order to obtain u q we need to analyze, again, the trace's argument in the expression 
for U q . We will follow the path already pursued above, but without integrating the q = 1-like 
term over frequencies. We obtain 

Tr (ZMH) = KAnJ* J duu n _ /Bto , (46) 

o e 

that, by recourse to the identity 

rrs;= £(«"*')• (47) 

allows us to cast U q in the fashion 



n = ^-i r [V(l q)\ hA 7 d n y yjl 7 dt ± — ■ , ( 4 f 

q ?ir I ^ ^ m / ft < mAt^ — ^ nm 

Z7T 5 s= 0m=0" L _ J m (l + Zt) 1 -i 

With Eq. ( p2[ ) the integral above is easily evaluated, and we obtain 
According to Eq. (p9|), the energy density defined by Eq. (£5^) is 

oo 

E 5 m ^ m r^ 1 [i/(i-g) + ^] 

= ft^a;"^ , (50) 

E 5 m r- 1 [l/(l-g) + nm] 

m=0 

where -B m is given by Eq. (|28|) while S" m reads 
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S m = J2[l-(l-q)P q huJs]^ +nm 



(51) 



where 



and 



q)(3 q huj_ 







(52) 



(53) 



l + (l-q)/3U q 

Notice that a new cut-off condition has been introduced, namely, 1 — (1 — q)/3 q hujs > 0, 
that transforms the original series into a finite sum. 

Since one knows U q from Eq. (SS), an expression for u q is easily obtained. Notice that it 
is not self referential. The shape of the resulting curves is similar to those arising from the 
traditional treatment, even for T values where one detects q- violations to Plank's law. Fig. 
f| depicts u q vs to for different q values. The maximum's values do not coincide with the 
ones of the canonical Gibb's approach when q differs from unity. In Fig. ^ (a) it becomes 
apparent that results for q = 1 and for q = 0.98 greatly differ. The difference seems to 
become less important between results for q = 0.9 and those for q = 0, as can be seen in 
(b). Another point to be stressed is that the maximum's position will also depend on q, not 
just on T, as in the traditional Wien's law. 

We have seen from Fig. |4] that Eq. (|50"D yields results that quite resemble the ones given 
by Planck's law. It is then reasonable to look for a perturbative expansion in 1 — q. The 
series over the m factor in Eq. fl5"U| ) rapidly converges. Only the first terms are important. 
In the limit q — > 1, the exponent g/(l — q) + nm ~ q/(l — q) in Eq. (|5l|) , and the finite sum 
becomes a series expansion. We have 



u q hA n uj n S, 



(54) 



where 



S='£[l-(l-q)f3 q hwsf: 



~I> X P [-q/3 q tuj(l + s)] 



(55) 
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a power series in s of guaranteed convergence, i.e., 



S ~ e q/3 q hu, Z 1 ' ( 56 ) 



so that replacement into Eq. fl54|) yields 



a first order correction to the Planck law. The classical result is attained for q — > 1. Eq. 
(|57|) provides one, then, with an approximate energy density. 

Following the standard text-book treatment p6f , we can define the particle-density as 

n q = (TiA n u n )~ l u qi (58) 

which, from Eq. ([57]) can be written, to first order in 1 — q as 

n « = (59) 

with e = hu. 

We have encountered an alternative expression for the boson particle-density of the black- 
body radiation, to be compared to the ones that result from either the FA or OLM-FA (see 
Apendix |A|). The correct value for the q — > 1 limit is obtained. 

Consider now the specific heat curves of Sect. [IV A| . From Eq. (B31) we see that, to first 
order in 1 — q we have 

[3 g ^f3{l-(l-q)f3U q ], (60) 

that, replaced into Eq. (|59|) and using the fact that q = 1 — (1 — q) leads to (keeping only 
terms of order 1 — q) 

where 

^ q = (I - q)e(l + (3U q ) (62) 
15 



plays the role of a fictitious g-chemical potential, such that //* — > for q — > 1. Of course, 
the "true" physical chemical potential vanishes identically (see Sect. pip . A sort of fictitious 
g-Bose-Einstein condensation effect might seem to be implied by the presence of this pseudo- 
chemical potential, a point that deserves further careful exploration, to be addressed in a 
future work. 



B. Wien law 

A nice result of the traditional thermostatistics is the linear relation that exists between 
the frequency u* that maximizes the energy density spectrum, on the one hand, and the 
temperature on the other one 

uj* oc T. (63) 

This relation is known as Wien's law. 

According to the results depicted in Fig. f| (see Sect. |V A] ), deviations of the these 
maximum values from what we expect from them according to the above relation will depend 
on q. If we perform the pertinent derivatives in Eq. (|50|) we obtain 



0, (64) 



£ BJT 1 [1/(1 -q) + nm] (-S m + ^ 

where S m depends in a non trivial manner on uj due to the cut-off condition (see Eq. (|TJ): 
the cut-off condition "cuts" the series at an s q that depends on to). 

Numerical calculations of Eq. fl6~4|) are depicted in Fig. [5]. It can be seen that the 
linear dependence is respected over temperatures above 1Q~ 2 K and below 10 _4 K. There 
is a "violation range" the T-zone (10~ 4 < T < 10~ 2 ). This implies lower values of T than 
those that violate Stefan Boltzmann's law. Departures from Wien's law are observed even 
for q— values as close to unity as q = 0.98. Using now Eq. ( p?|) as a starting point we also get 
first order corrections in 1 — q to the Wien law. We can locate the maximum with respect 
to uj of Eq. ( p?| ) with the auxiliary definition 
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ujq(3h 



(65) 



x = 



1-(1-<Z)/HV 



so that we immediately find 




(66) 



whose solution is a constant, 6, for each fixed value of n. For n = 3 we find, for instance, 
b = 2.82. The maximum of u q {ui) for different T's is located at distinct frequencies according 
to 



The OLM formalism can be applied to the analysis of the cosmic microwave radiation in 
order to search for putative deviations from Planck's law. The best source of data on this 
topic comes from the COBE satellite (Far-infrared Absolute Spectrophotometer (FIRAS)). 
The ensuing brightness, when compared to a Planck spectrum at the temperature T = 
2.72584 ± 0.0005A" shows significant deviations from the expected form (the x 2 °f the fit to 
a Planck spectrum is rather large). 

FIRAS measures the differential spectrum between the cosmic background and an in- 
ternal reference adjusted to approximately 2.7A" so as to i) avoid problems with imperfect 
emissivity and ii) limit the dynamic range of the instrument. In addition, an external black- 
body is used to calibrate (according to Gibbs' statistics) the gain of the instrument with a 
temperature in the range of IK up to 25A". 

This is the kind of data that will be of interest in order to verify the predictions of the 
OLM formalism. If we assume the cosmic background does not obey the usual Planck's 
law, but the OLM one, we can check the concomitant differences and contrast the ensuing 
values with those obtained by FIRAS. The results are depicted in Fig. || We plot the 




q h ' 



(67) 



a first order correction to Wien's law. 



C. Microwave Cosmic Background Radiation 
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best fits for four different values of q. The range of temperatures that allow for the best 
fitting is also given. In all the graphs it is apparent that deviations are well predicted if we 
assume a statistics with q ^ 1 for the cosmic background. As q deviates from unity higher 
temperatures are predicted, in agreement with the results displayed in Fig. ^. We see there 
that the spectrum maximum's value becomes smaller for a given, fixed temperature when q 
deviates from unity. 

The experimental data depicted in Fig. |] seems to imply that the statistics ruling the 
cosmic microwave radiation is different from that of the conventional black-body instance. 
This constitutes evidence for the existence of alternative statistics, and tells us that the 
statistics behind these phenomena seems to possess an additional, g-degree of freedom. 

VI. CONCLUSIONS 

Within the g-Thermostatistics framework we perform an exact analysis of the n- 
dimensional black-body radiation process. We employ to such effect both Tsallis' and Renyi's 
entropies, within the range < q < 1. The new theoretical ingredient here is the so-called 
OLM approach to non-extensive thermostatistics [|TjJ. 

We develop a q generalization of several laws: Stefan Boltzmann's , Planck's, and Wien's. 
We find the the conventional behavior still obtains for temperatures above IK, although 
there is a q dependence in the appropriate proportionally constants. We recover the tradi- 
tional relationship between radiation pressure and internal energy using the OLM formalism. 

We apply the formalism to experimental data on the cosmic microwaves' background and 
reproduce it with acceptable accuracy even for g-values that appreciably differ from unity. 
The larger |1 — q\ is, the higher the predicted equilibrium cosmic temperature. 
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1. Un-normalized Factorization Approximation (FA) 



In order to facilitate the reader's task we review now the factorization approach (FA) 
treatment, due to Biiyiikkilic et al., who tackled in Ref. |23[] the quantum ideal gas problem 
in a grand canonical scenario using the CT-Tsallis formalism. Tsallis' entropy is thereby 
maximized subject to the un-normalized constraints 

Trp q H = U q (Al) 
Trp q N = N q . (A2) 

The pertinent partition function is 

^ r = E[i-(i-ff)E^l 1/(1 " 9) . (A3) 

i 3 

where n i3 - are the occupation numbers of the level j (with single particle energy 6j) for a 
given i-microscopic configuration and 

x' 3 =f5 CT {e 3 -^ T ), (A4) 

with (3 CT the inverse temperature and p CT the chemical potential. 



According to [23|, "there is no restriction on the summation and thus one an factorise 
a product of factors, one for each one-particle state, since particle are regarded statistically 
independent" . Therefore (the essential point of the FA approach), the partition function 
Zq T given by Eq. (p3|) can be factorized. Each factor corresponds to a single particle state 



oo to 



? T -nE[i-(i-M 1/(H) ' (A5) 

j=0 i=0 



The average occupation number {n,j) FA of the state j becomes then [[23 

{U])fA ~ [l-(l-g)^]V(l-«) T l' (A6) 

which is a non-extensive, FA-Tsallis generalization of the Bose-Einstein (Fermi-Dirac) dis- 
tribution. 
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2. The normalized OLM treatment 



An alternative treatment to the factorization approach using the TMP formalism yields 
self-referential probabilities. The OLM procedure, instead, is not afflicted by such a problem. 
One again maximizes Tsallis' generalized entropy given by Eq. (|TJ) subject to the diagonal 
normalized constraints for the Grand Canonical Ensemble, which results in a partition func- 
tion of the form 



!-(!-?) E 



1/(1-9) 



(A7) 



where e x q = [1 + (1 — q)x] 1 ^ 1 q ^ and, instead of x'j (Eq. 



)), we have 



(A8) 



with 



x q = p(U g - (XNg). 



(A9) 



Note that Xj — > x'j for q — > 1. The conventional limit is the same for both treatments, if 
one remembers that the un-normalized methodology uses a partition function Z q - !T , while the 
normalized treatment deals with a bar partition function Z q (see Eq. (|TT|)). The concomitant 
q = 1-limits are related according to Z% = e~^ v ~ tlN ' Z\. 

Proceeding once again as in Ref. PB| , i.e., neglecting correlations between particles and 
regarding states of different particles as statistically independent, the partition function Z q 
can be factorized. Each factor corresponds to a single particle state, 



^-<'nE[i-(i-^»A'] 1/M 



(A10) 



3 « 



The present expression for Z q is formally identical to that encountered in Ref. []23| (see Eq. 
J)), save for a multiplicative factor, which does not contribute to the ensuing microscopic 
probabilities. Thus, according to Ref. the average value (ti^nfa is now 

1 



( n i) 



NFA 



[1 - (1 - q)x j ] 1 /( 1 -i) T 1 
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(All) 



with Xj given by Eq. 

We introduce here the normalized factorization approach (NFA) black-body's treatment. 
The mean occupation number is given by Eq. (|A11| ) with fi = (the chemical potential 
vanishes because we are working with photons). In the continuum limit e(uj) = Tiu (u is the 
frequency), and the particle-density will be given by 

{n{u})) ^ w [i-(i- g yHF/(H)-r (A12) 

with 

x M = i + a-w (A13) 

In the three dimensional case the mean occupation number is connected with the energy- 
density according to 

u q (u) = ^-y(n(u)) g . (A14) 



Keeping only first order terms in 1 — q in Eq. (|A14 ), integrating over frequencies, and 
re-arranging terms, the internal energy reads 



AT. \ , 1 ( 1 - <l ]0 

l + *(l-q)o7« 



jjNFA = aT 4 _ - _±± (A15) 



where a = (7r 2 A; 4 V)/(15/i 3 c 3 ) (with c, the light's speed) is in closed connection with the 
Stefan-Boltzmann constant (ccr) / (4V) . 
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FIGURES 
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FIG. 2. Generalized Stefan constant a q as a function the index q. It is shown \n(a q /a) vs q 
with a = 5.67051 • 10" 8 Watt/(i^ 4 Meter 2 ) the Stefan constant for g = 1. See inline text for details. 
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FIG. 3. Specific heat, C q , of a three-dimensional box of volume lm 3 , as a function of Tempera- 
ture T, for different values of the non-extensivity parameter q. The plots show the arising results 
for (a) q = 0.98, (b) q = 0.8, (c) q = 0, (b) q = 1. 
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FIG. 4. Energy density n 9 as a function of the frequency u>. All the curves were obtained for a 
temperature T = O.li^, and V = lm 3 . Inset (a) shows the results for the limit q — > 1, whereas in 
the inset (b) the results for obtained otherwise are depicted. 
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FIG. 5. Frequency uj* for which the energy density reaches its maximum, as a function of 
temperature T, for different values of q. As in the previous figures, n = 3. 
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FIG. 6. We show the differential radiation between for the FIRAS spectrographer onboard 
COBE satellite. Superimpossed with the experimental data the best fits for four different q values 
are also shown. 
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TABLES 



Formalism 


Result 




Stefan-Boltzmann' law 


Ui = aT A 


a = (tt 2 A; 4 \/)/(15^ 3 c 3 ) 


C-T Solution 


IJCT = 3£ 3 kTT l -o 


[2-,- 


E f r-^[2Ef+3m+3] 

m = 


„ _ 4r(3)C(4)V / v^fcT \ 3 


U-5. 




E ^r-i[|5f+3 m ] 

m = 




« - r(3/2) Uc(i-?)J 


OLM Solution 


CO 

, J2 Bmr- 1 [ T l-+3m+4l 

ttOLM,R Qt U rr[l+{l-q)PU q ] m=0 


d _ [l+(l-g)/3C/ a r m 


u q — °?3«' J (1-a) 4 

E B m F-^ + 3m] 




FA Aproximation 
OLM-FA Aproximation 


U( FA > « <r[l - (1 - q)6}T 4 

jtNFA ^ _ T 4 
u g ~ a 1+4(1- ? )ctT3 /fe J 


61 = C(5)r(6)/(2C(4)r(4)) 


OLM Aproximation 


ttOLM 1 7^4 

~ u i-4(i- 9 ) CT T 3 /fe 





TABLE I. The results of different generalizations of Stefan-Boltzmann's law. 
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